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Abstract 


This  work  presents  a  circuit  model  for  calculating  the  total  energy  dissipated  into  neutral  species 
for  pulsed  direct  current  (DC)  dielectric  barrier  discharge  (DBD)  plasmas.  Based  on 
experimental  observations,  it  is  assumed  that  nanosecond  pulsed  DBDs,  which  have  been 
proposed  for  aerodynamic  flow  control,  can  be  approximated  by  the  two  independent  regions  of 
a  homogeneous  electric  field.  An  equivalent  circuit  model  is  developed  for  the  homogeneous 
region  near  the  exposed  electrode,  i.e.,  the  “hot  spot,”  based  on  a  combination  of  a  resistor, 
capacitors,  and  a  zener  diode.  Instead  of  fitting  the  resistance  to  an  experimental  data  set,  a 
formula  is  established  for  approximating  the  resistance  by  modeling  a  plasma  as  a  conductor 
with  DC  voltage  applied  to  it.  Various  assumptions  are  then  applied  to  the  governing  Boltzmann 
kinetic  energy  equation  to  approximate  electrical  conductivity  values  for  weakly  ionized 
plasmas.  The  model  is  compared  with  experimental  data  sets  of  the  total  power  dissipated  by  a 
plasma  to  validate  its  accuracy. 


154 


Acknowledgments 


I  wish  to  thank  my  mentor,  Dr.  Bryan  Glaz,  for  the  guidance  and  support  he  provided  throughout 
the  period  of  research.  I  would  also  like  to  acknowledge  the  support  provided  by  Dr.  Surya 
Dinavahi  and  Tomas  Houba  regarding  the  direction  of  research  as  well  as  the  many  enlightening 
discussions  shared  with  them. 


155 


Student  Bio 


I  am  currently  an  undergraduate  student  attending  the  University  of  Florida  (UF),  majoring  in 
both  nuclear  engineering  and  physics.  At  UF,  I  am  a  member  of  the  Applied  Physics  Research 
Group  headed  by  Prof.  Subrata  Roy,  where  I  conduct  research  on  plasma  physics.  In  May  of 
201 1, 1  completed  his  second  year  of  studies  at  UF  and  was  lucky  enough  to  have  an  opportunity 
to  work  with  plasma  actuators  on  a  project  of  interest  for  the  U.S.  Army.  Upon  completion  of  my 
undergraduate  education,  I  hope  to  pursue  the  subject  of  applied  plasma  physics  in  graduate 
school  in  order  to  obtain  a  doctoral  degree. 


156 


1.  Introduction/Background 


The  need  for  improved  control  over  aerodynamic  flow  separation  has  increased  interest  in  the 
potential  use  of  plasma  actuators.  The  inherent  advantages  of  plasma  actuator  flow  control 
devices  include  minimal  size  and  weight,  lack  of  moving  parts,  cost  efficiency,  and  low  drag 
penalties.  However,  the  actuators  that  affect  the  flow  via  directed  momentum  transfer  are  not 
effective  at  Mach  numbers  associated  with  most  subsonic  aircraft  applications.  Recently, 
Roupassov  et  al.  ( 1 )  demonstrated  that  pulsed  plasma  actuators,  in  which  energy  imparted  to  the 
flow  appears  to  effectively  control  flow  separation,  seem  to  be  suitable  at  Mach  numbers  (M  ~ 
0.3)  beyond  the  capabilities  of  the  current  plasma-induced  momentum-based  approaches. 

Given  the  fundamental  differences  between  the  novel  pulsed  discharge  approach  and  the  more 
conventional  momentum-based  approaches,  there  is  a  need  to  develop  an  effective  and  efficient 
model  for  the  energy  delivered  to  the  flow  by  the  plasma.  Once  calculated,  that  value  can  be 
input  to  a  computational  fluid  dynamics  solver  as  an  energy  source  term  resulting  in  a  coupled 
fluid/plasma  dynamics  model.  Multiphysics  models  of  this  type  are  required  in  order  to  study 
detailed  flow  characteristics.  However,  detailed  numerical  calculations  that  are  primarily  used 
for  isolated  plasma  simulations  are  not  suitable  for  a  variety  of  coupled  fluid/plasma  dynamic 
studies  because  of  their  excessive  computational  expense.  To  address  this  issue,  efficient  circuit 
element  models  have  been  introduced  to  approximate  the  complex  dynamic  processes  within 
plasmas.  However,  models  such  as  those  by  Orlov  et  al.  (2)  rely  on  empirical  constants  that  are 
tuned  to  experiments  which  are  not  applicable  to  nanosecond  pulsed  discharges.  To  date,  an 
approximate  model  of  nanosecond  pulsed  plasma  actuators  has  not  been  developed.  This  paper 
deals  primarily  with  establishing  a  flexible  model  that  could  be  implemented  as  an 
approximation  for  the  energy  dissipated  within  a  plasma  for  any  pulsed  direct  current  (DC) 
dielectric  barrier  discharge  (DBD)  configuration.  Among  the  other  fundamental  goals  in  this 
paper  is  to  establish  a  deeper  understanding  of  the  background  processes  that  occur  within 
plasma  and  incorporate  that  knowledge  into  the  model. 


2.  Theoretical  Development 


One  of  the  primary  assumptions  in  creating  this  model  is  that  a  small  volume  near  the  exposed 
electrode  can  be  idealized  as  consisting  of  a  homogeneous  energy  density.  The  “hot  spot,” 
observed  by  Roupassov  et  al.  (1 ),  plays  an  important  role  in  the  power  dissipated  into  the  air  and 
thus  the  resulting  fluid  dynamics.  In  this  study,  only  the  spatially  homogeneous  hot-spot  region 
was  considered.  Therefore,  the  model  presented  in  this  paper  consists  of  a  single  network 
containing  a  resistor,  capacitors,  and  a  diode. 
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As  shown  in  figure  1 ,  circuit  elements  that  were  used  to  model  the  plasma  include  an  air 
capacitor  Ca,  a  dielectric  capacitor  O,  a  resistor  Rj,  and  a  zener  diode  Df.  The  air  capacitor 
represents  the  capacitance  between  the  dielectric  surface  and  the  exposed  electrodes.  The 
dielectric  capacitor  represents  the  capacitance  between  the  dielectric  surface  and  insulated 
electrodes  and  is  proportional  to  the  thickness  of  the  dielectric  layer.  Thus  the  dielectric  layer  in 
the  form  of  both  its  thickness  and  the  value  of  its  dielectric  constant  plays  an  important  role  in 
determining  the  effectiveness  of  the  plasma  actuator.  Finally,  the  zener  diode,  introduced  by 
Orlov  et  al.  (2),  is  used  in  the  model  to  enforce  an  energy  threshold  value  below  which  plasma 
will  not  form. 


Figure  1.  Electric  circuit  model  of  a  dielectric  aerodynamic  plasma  actuator. 

Since  a  uniform  charge  distribution  along  the  top  of  the  dielectric  is  assumed,  the  typical 
asymmetric  plasma  actuator  geometry  featured  in  figure  1  can  be  simplified  to  a  homogeneous 
symmetric  region.  This  assumption  results  in  a  one-dimensional  (1-D)  model,  and  any  variation 
along  the  horizontal  direction  depicted  in  figure  2  is  neglected. 


Figure  2.  Region  of  homogeneous 
potential,  i.e.,  “hot  spot.” 

As  displayed  in  figure  1,  the  lumped  element  circuit  is  a  function  of  the  two  capacitance  values, 
Ca  and  Cd-  In  this  model,  the  air  is  treated  as  both  a  conductor  to  generate  a  physical  relationship 
for  the  resistance  and  a  parallel  plate  capacitor  to  generate  Ca.  An  advantage  of  modeling  the 
plasma  as  a  conductor  in  addition  to  a  parallel  plate  capacitor  is  that  it  generates  a  physical 
relationship  for  the  resistance,  Rf — a  value  that  is  traditionally  empirically  determined.  The  air 
gap  capacitor  can  be  modeled  as  (3) 
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where  Aa  is  the  cross-sectional  area  of  the  air  and  ha  is  the  approximate  height  of  the  plasma 
region  of  interest.  As  displayed  in  figure  3,  Aa  is  the  product  of  the  spanwise  length  of  the 
actuator  za,  and  la  is  the  chordwise  distance  from  the  exposed  electrode  to  the  end  of  the 
dielectric  region. 


Exposed  Electrode 


Air  Capacitance  Element 


Dielectric 


•VdAVA  iy'f-U«o!r 


e_ 


Figure  3.  Sketch  of  the  air  capacitive  element. 


The  capacitive  element  corresponding  to  the  dielectric  can  be  modeled  as  (3) 

r  _  e0  edAd 

Cd~~> 


(2) 


where  Ad  is  the  cross-sectional  area  of  the  dielectric  capacitive  element  and  hd  is  the  height  of  the 
dielectric  barrier  layer.  As  displayed  in  figure  4,  Ad  is  the  product  of  the  spanwise  length  of  the 
actuator  za  and  dd  is  the  width  of  the  dielectric  region.  Treating  the  plasma  as  a  conductor,  the 
resistance  for  DC  voltage  is  proportional  to  the  electrical  conductivity,  op,  as  well  as  Aa,  and  ha 
and  can  be  given  as  ( 4 ) 


hq 

<JpAa 


(3) 


Figure  4.  Sketch  of  the  dielectric  capacitive  element. 

Starting  from  Kirchhoff  s  laws  (3),  the  governing  differential  equation  for  the  voltage  drop 
experienced  by  the  air  gap  AC  is  given  by 

dAV{£)  _  _  dVgppjt )  /  ca  _  \  _  A v-(t) 

dt  dt  \Ca+Cd  )  «/(t)(Ca+Cd) 
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and 


(nf\E\>mcrit 

(o  if  \E\  <  \E\crit 


where  Vapp  is  the  applied  voltage  and  k  is  the  contribution  from  the  zener  diode.  If  the  electric 
field  magnitude,  given  as 


(6) 


is  greater  than  some  threshold  value,  then  k  takes  on  a  value  of  one,  otherwise  it  is  zero  to  signify 
that  plasma  has  not  formed.  For  nanosecond  high-voltage  pulses,  the  applied  electric  field  will 
virtually  always  be  over  the  critical  breakdown  field  (7 ),  thus  k  can  be  set  equal  to  one  for  the 
case  considered  in  this  study. 


The  energy  imparted  to  neutral  species  from  the  plasma  can  be  calculated  from  knowledge  of  the 
voltage  drop  over  the  air  gap  by 


E(t)  — 

V  J  J0  Rf(t ) 


(7) 


To  effectively  calculate  the  resistance  governed  by  equation  3,  an  expression  must  first  be 
developed  for  the  electrical  conductivity  of  the  plasma.  This  value  is  one  that  traditionally 
requires  a  numerical  approach.  To  simplify  the  problem  to  a  point  where  an  analytic  fonnulation 
can  be  used,  numerous  simplifying  assumptions  were  used  and  are  described  in  the  following 
paragraphs. 


For  any  plasma,  the  resulting  electric  current  is  composed  of  two  primary  terms:  the  current  from 
electrons  and  that  from  ions.  As  the  drift  velocity,  we,  which  represents  the  velocity  induced  by 
an  electric  field,  is  significantly  higher  in  a  nonequilibrium  plasma  for  electrons  compared  to 
ions,  the  current  density  can  be  approximated  as  only  the  portion  from  electrons  as  long  as  the 
number  densities,  Ne  and  TV*,  are  approximately  the  same  ( 4 ).  When  a  fonn  of  the  generalized 
Ohm’s  law  is  used,  the  current  density  vector  /  and  plasma  conductivity,  respectively,  can  be 
written  as 


/  w  —eNewe  =  <rpE, 
and 

op  =  e(Nefie  +  NiHi), 


(8) 

(9) 


where  /q  and  jue  represent  the  ion  and  electron  mobilities,  respectively.  Much  like  equation  8,  the 
electrical  conductivity  relation  can  be  simplified  using  the  concept  of  quasineutrality,  which  is 
defined  as  having  approximate  equal  number  densities  of  both  ions  and  electrons.  Thus,  as  iie  is 
typically  approximately  three  orders  of  magnitude  larger  than  //  ,  for  any  plasma  that  is  not  in 
thermal  equilibrium,  if  Ne  is  at  least  the  same  order  of  magnitude  as  Ni,  it  is  a  good  assumption  to 
approximate  the  electrical  conductivity  as  only  coming  from  electrons  ( 4 ).  Quasineutrality  itself 
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is  a  typical  assumption  that  is  valid  as  long  as  the  plasma  being  modeled  is  far  away  from  the 
cathode  to  avoid  the  boundary  layer  in  plasma  physics  called  the  sheath. 

Since  a  pulsed  DC  voltage  is  assumed,  the  activation  of  the  external  electric  field  will  follow  the 
voltage  waveform  as  a  step  function.  Thus  two  expressions  will  be  required  for  the  op,  where  the 
first  is  valid  for  the  period  when  a  constant  external  electric  field  is  applied  over  the  peak  of  the 
voltage  pulse,  as  shown  in  figure  5,  from  13-40  ns,  and  the  second  when  the  external  voltage  is 
zero.  For  the  portions  of  the  voltage  wavefomi  that  the  voltage  is  zero,  the  power  is  also  zero 
according  to  Ohm’s  law,  and  thus  the  conductivity  during  this  time  is  of  no  importance. 


Figure  5.  Input  voltage,  Vapp  vs.  time. 


Given  an  assumption  of  a  constant  electric  field  during  the  peak,  an  analytic  formulation  can  be 
generated  for  the  electrical  conductivity  with  knowledge  of  we  (5)  by  starting  with  the 
Boltzmann  equation  for  electrons  written  as 


d±  +  v.d-L  +  L.d-L=V' 

dt  dr  m  dr  dt’ 


(10) 


where  df/dt  is  typically  called  the  collision  integral  and  refers  to  the  variation  in  the  number  of 
particles  within  some  phase  space  control  volume.  In  order  to  approximate  this,  a  relaxation 
time,  r,  is  introduced,  which  is  defined  as  the  time  taken  for  the  system  to  be  reduced  to  an 
equilibrium  distribution  function.  This  approximation  can  be  given  as  t  ~  (Nacreav)"  ,  where  oea 
is  the  collision  cross  section,  v  is  the  average  collision  velocity,  and  Na  is  the  number  density  of 
atoms  (5).  According  to  Smirnov  (5),  an  analytic  formulation  is  obtained  by  assuming  that  the 
force  term  can  be  given  as  the  Lorentz  force, 

F  —  —  eEexp(—  icot)  —  ^  (v  x  B),  (11) 


where  a)  is  the  frequency  the  electric  field,  B  is  the  magnetic  field  and  assumed  to  be  directed 
along  the  z-axis  (vertical),  and  E  is  assumed  to  belong  only  in  the  x(horizontal)-z  plane.  Taking 
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the  tau  approximation  of  electron  kinetic  energy  equation  and  integrating  over  the  electron 
velocity,  we  can  then  write  the  equation  of  motion  for  an  electron  as 

me  ~  +  me  ^  exp(-icot) — (v  x  B).  (12) 


Substituting  the  cyclotron  frequency  o)H  —  eB/mec,  ax  —  —eEx/me  and  az  —  —eEz/me  ,  the 
following  scalar  equations  representing  the  different  components  of  equation  12  can  be 
produced: 


-jj-  +  —  =  ax  exp  (—nut)  —  ooHwy, 

dWy  Wy 

—  +  T  =  0)H0)X, 
and 


(13) 

(14) 


dwz  wz  ,  .  A 

— —  +  —  =  az  exp(— loot). 


(15) 


Equations  13-15  can  be  solved  to  obtain  three  component  expressions  for  the  electron  drift 
velocity.  If  a  constant  electric  field  and  a  negligible  magnetic  field  are  assumed,  a>  and  odH  can 
be  set  equal  to  zero  and  the  resulting  vector  equation  can  be  obtained: 


eEx 


(16) 


When  equation  16  is  solved  in  conjunction  with  equation  8  for  erp,  the  following  expression 
independent  of  the  applied  electric  field  can  be  obtained: 


We(t)e2 

meNaaeav' 


(17) 


In  equation  17,  oea  is  a  function  of  the  electron  temperature,  Te,  and  can  be  obtained  for  various 
molecules  found  in  air  (<5).  As  shown  in  equation  17,  the  inputs  required  for  this  approximation 
include  Ne  ,  the  number  density  of  electrons,  Na,  and  Te.  Among  these  values,  Na  is  assumed  to 
be  constant  in  time  as  the  number  density  of  atoms  is  significantly  higher  than  that  of  free 
electrons.  Te  is  solved  for  quantity  from  the  energy  transport  equation  and  changes  primarily 
with  reduced  electric  field  strength,  E/Na.  Results  for  the  electron  temperature  are  generated 
from  BOLSIG+,  a  program  aimed  at  solving  for  the  transport  properties  of  electrons. 


The  final  required  input,  the  electron  number  density  Ne,  is  traditionally  solved  from  the  drift 
diffusion  approximation  of  the  conservation  of  momentum  equation.  Starting  from  the  species 
continuity  equation  for  electrons  and  neglecting  diffusion,  an  equation  for  the  time  rate  of  change 
in  electron  density  can  be  obtained: 

^=a\WeE\-pn2,  (18) 


where  a  determines  the  ionization  frequency  and  /?  establishes  the  recombination  rate  of 
electrons.  As  equation  18  is  function  of  the  electric  field  experienced  by  a  plasma,  a  complete 
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system  of  equations  will  require  equations  4,  6,  17,  and  18  solved  in  a  coupled  manner.  Equation 
18  is  also  a  function  of  empirical  coefficients  a  and  /?  which  are  both  dependent  on  the  pressure, 
p,  of  the  system  and  can  be  given  by 

a  =  Ap  exp  (~jp  ),  (19) 

and 

p = c  (™y,  (20) 

where  A,  B,  and  C  are  tabulated  values  for  air  at  atmospheric  pressure.  According  to  Raizer  (7), 
appropriate  values  based  on  experimental  measurements  are:  38.1  (in-Torr)'1  (15  [cm-Torr]'1 )  for 
A,  927.1  V/(in-Torr)  (365  [V/(cm-Torr)])  for  B  and  1.2  x  10-8  in’/s  (2  x  1CT7)  cm3/s  forC. 


3.  Results  and  Discussion 


In  order  to  validate  the  accuracy  of  the  model  described  in  this  paper,  comparisons  with  data 
presented  in  Roupassov  et  al.  (7)  are  provided.  The  experimental  parameters  that  were  mentioned 
and  used  in  the  circuit  model  are  given  in  table  1 . 


Tablet.  Experimental  parameters. 


h-a 

0.016  in  (0.4  mm) 

hd 

0.01  in  (0.3  mm) 

ed 

2.7 

ea 

1 

An 

0.05  in'  (30  min2) 

Ad 

0.05  in'  (30  mm3) 

V 

50  kV 

T 

50  ns 

An  approximation  of  the  50-kV  applied  voltage  pulse  employed  by  Roupassov  et  al.  (7)  was 
considered  for  this  study  and  is  shown  in  figure  5. 

For  an  E /Na  value  of  4921  Td,  BOLSIG+  calculates  that  the  mean  electron  energy  should  be 
~1.18  x  10_2°  Btu  (78  eV).  Using  this  result,  we  calculated  a  peak  number  density  of  ~6. 5  x 
1017 1/ft3  (2.3  x  1019l/m3)  using  equation  18  coupled  with  equation  4  where  7V0,  the  number 
density  before  the  pulse,  was  set  equal  to  2.83  x  1013 1/ft3  (1015  1/m3).  The  height  of  the 
plasma,  ha,  was  determined  to  be  a  constant  0.02  in  (0.4  mm)  over  time  based  on  experimental 
observations  by  Roupassov  et  al.  (7).  Thus  using  the  calculated  values  for  Ne(t )  and  Te  as  well 
as  the  experimentally  observed  ha,  we  obtained  a  value  of  2.3  x  10~6  Btu  (2.6  mJ)  by  using 
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equation  7.  Gas  heating  results  from  Roupassov  et  al.  (7)  correspond  to  a  total  energy  of  4.2  mJ 
that  was  imparted  to  air  per  65  ns.  Therefore,  the  approximate  model  presented  in  this  study 
provides  an  answer  with  an  absolute  error  of  -38%.  The  accuracy  of  higher-fidelity  numerical 
models  based  on  discretization  of  the  drift  diffusion  equations  has  not  been  established. 


4.  Summary  and  Conclusions 


A  new  lumped  element  circuit  model  was  presented  that  is  valid  for  any  pulsed  DC  DBD  plasma. 
In  addition,  an  approximate  expression  was  fonnulated  to  calculate  the  resistance  value  for  the 
air  gap  as  a  function  of  the  conductivity  of  the  plasma.  Results  of  the  model  were  verified  against 
a  pulsed  DC  experiment  conducted  by  Roupassov  et  al.  (7),  and  an  order  of  magnitude 
agreement  was  obtained  for  the  energy  imparted  into  the  plasma  in  a  homogeneous  “hot  spot” 
region.  Future  work  will  focus  on  extending  the  circuit  model  to  account  for  the  plasma  region 
outside  of  the  “hot  spot”  as  well  as  determining  the  geometric  dimensions  of  the  plasma  regions 
as  functions  of  the  applied  voltage. 
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